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Abstract 

We construct stable axially symmetric models of elliptical galaxies. 
The particle density on phase space for these models depends mono- 
tonically on the particle energy and on the third component of the an- 
gular momentum. They are obtained as minimizers of suitably defined 
energy-Casimir functionals, and this implies their nonlinear stability. 
Since our analysis proceeds from a rigorous but purely mathematical 
point of view it should be interesting to determine if any of our mod- 
els match observational data in astrophysics. The main purpose of 
these notes is to initiate some exchange of information between the 
astrophysics and the mathematics communities. 



1 Introduction 

Consider a large ensemble of mass points (stars) which interact only by the 
gravitational field which they create collectively. Such a collisionless, self- 
gravitating gas is used to model galaxies or globular clusters, cf. p, ^] and 
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the references there. The time evolution of the density / = f(t,x,v) > of the 
stars in phase space is governed by the following nonlinear system of partial 
differential equations: 

d t f+vd x f-d x U-d v f = 0, 
AU = Att P , lim U{t,x) = 0, 

\x\^oo 

p(t,x) — J f(t,x,v)dv. 

Here tGlR denotes time, x,v GlR 3 denote position and velocity respectively, 
p is the spatial mass density of the stars, and U the gravitational potential 
which the stars induce collectively. Collisional or relativistic effects are ne- 
glected. This system, which Sir J. H. Jeans introduced at the beginning of 
the last century for the above modeling purposes, is sometimes in the astro- 
physics literature referred to as the collisionless Boltzmann-Poisson system. 
Following the mathematics convention we call it the Vlasov-Poisson system. 

In the present paper we are interested in the steady states of this system. 
Besides their very existence, questions of interest are the possible shapes of 
such steady states, in particular their symmetry type, and whether or not 
they are dynamically stable. We approach these questions from a mathemat- 
ics point of view, in particular, only such information as can be extracted 
from the system stated above is to enter our arguments. Our "ideal" reader 
is an astrophysicist interested in these same questions. We find it deplorable 
that there is little communication between astrophysicists and mathemati- 
cians investigating these problems, and the present paper is an attempt to 
change this. Axially symmetric steady states of the the Vlasov-Poisson sys- 
tem were obtained in as perturbations of spherically symmetric ones via 
the implicit function theorem. From an astrophysics point of view axisym- 



metric models of elliptical galaxies have been investigated in [ I8| , |2T| , p8fl , but 
as to their stability no rigorous results are known to us. In the present paper 
we follow the variational approach developed in 0, §, g ^ [FT], JZQ We 



obtain axially symmetric steady states with vanishing and with non- vanishing 
velocity field as minimizers of appropriately defined energy-Casimir function- 
als. These steady states are no longer restricted to small perturbations of 
spherically symmetric ones, and, most importantly, they are nonlinearly sta- 
ble. The question we would like to raise is: Do our mathematical construc- 
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tions yield suitable models for real- world galaxies, or, if not, how can the 
mathematical approach be modified to do so? 

If Uq = U q {x) is a stationary potential then any function f = (p(E) of the 
particle energy 

E = E(x,v) = ~\v\ 2 + U (x) (i.i) 

solves the Vlasov equation with the potential Uq, since E is a conserved 
quantity along particle trajectories. So if U is the potential induced by f 
then this is a steady state of the Vlasov-Poisson system — at this point it 
must be emphasized that to obtain a self-consistent model by this approach, 
i. e., to make sure that Uq is indeed the potential induced by fo, is mathe- 
matically non-trivial. Any steady state obtained by this ansatz is necessarily 
spherically symmetric, as follows from a result of Gidas, Ni & Nirenberg |J. 
In the present paper we make fo depend on an additional invariant of the 
particle trajectories. If Uq is axially symmetric, i. e., invariant under all ro- 
tations about, say, the 23-axis, then the corresponding component of angular 
momentum, 

P = P(x,v) = xiv 2 — x 2 Vi (1.2) 

is such an additional invariant. To obtain steady states with a number density 
fo which depends on E and P we minimize the energy-Casimir functional 

\JJ \v\ 2 f(x,v)dvdx-^J J dxdy+ J I Q(f(x,v),P(x,v))dvdx 

subject to the constraint that / is axially symmetric and has prescribed mass 
M>0: 

f(x,v)dvdx = M. 



IS- 



The corresponding Euler-Lagrange equation then shows that any minimizer 
is a steady state of the desired form. In order to conclude that this steady 
state is nonlinearly stable it is essential that the energy-Casimir functional 
is conserved along solutions. This is always true for the sum of the first two 
terms which are the kinetic and the potential energy of the system. However, 
the third term is conserved only if P is constant along particle orbits which 
is true in general only if the solution is axially symmetric. Thus, the stabil- 
ity result which we derive is restricted to axially symmetric perturbations of 
the steady state. The fact that the third term in our functional is not con- 
served along general solutions makes the terminology "Casimir functional" 
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questionable, but we stick to it for lack of a better alternative. The situation 
is similar to the one of anisotropic, spherically symmetric steady states with 
f depending on the particle energy and the modulus of the angular mo- 
mentum, \x X v\, where stability against spherically symmetric perturbations 
holds, cf. [flOfl , but the general case is open. For isotropic steady states where 
fo depends only on the particle energy no such restriction is necessary, cf. 
0,11- 

We want to stress the fact that in this article every conclusion is rigorously 
derived in the mathematical sense. Our method stems from our long term 
research project in which both the existence and the stability of galaxy con- 
figurations has been mathematically investigated. Since the Vlasov-Poisson 
system can be viewed as an infinite dimensional Hamiltonian system, even 
the celebrated Antonov stability theorems for polytropes are, strictly speak- 
ing, only valid in the linearized sense. This is like the case of center equilibria 
in a 2 x 2 Hamiltonian system: in general no dynamical stability assertion 
can be deduced from a linearized analysis, and some Liapunov type meth- 
ods are needed. Recently, we were able to verify the validity of Antonov's 
theorems in a dynamical, nonlinear sense [pp . Our variational method also 



yields stability of certain Camm models [1C . 

The paper proceeds as follows: In the next section we formulate the basic 
set up and state our main results. In order to put these into perspective 
we include a short discussion of stability concepts from a mathematics point 
of view. In the third section we discuss various examples and some further 
properties of the steady states which we obtain. For interested readers we 
collect the details of our mathematical proofs in an appendix. 

We conclude the introduction with some references to the mathematical 
literature. In addition to our work mentioned above the stability of spher- 
ically symmetric steady states of the Vlasov-Poisson system in the present 
stellar dynamics case is also investigated in ^ |2{J 0. Global classical so- 
lutions to the initial value problem for the Vlasov-Poisson system were first 
established in W%, cf. also pfl and ||20||. For the plasma physics case where 



the sign in the Poisson equation is reversed, the stability problem is better 
understood, and we refer to ML O, [13, pL9| . A rather general condition which 



guarantees finite mass and compact support of steady states, but not their 



stability, is established in p5 
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2 Main results 

Let us first fix some notation. For a measurable function f = f(x,v) >0 we 
define its induced spatial density 

Pf(x):= J f(x,v)dv, xGlR 3 , 

and its induced gravitational potential 

u f'-=-Pf*yr 

Next we define 



E ki n{f) := \JJ \v\ 2 f{x,v)dvdx, 



EUf) ■= ~ J \VU f (x)fdx = - 1 -J J P^^ dxdy, 



C(f) := J J Q(f(x,v),P(x,v))dvdx 
where Q is a given function satisfying certain assumptions specified below 



and P is defined in (|1.2|) . We will minimize the energy- Casimir functional 

'He '■ = Ekin + E pot + C 
under a mass constraint, i. e., over the set 

^m:={/G^(IR 6 )|/>0, J Jf(x,v)dvdx = M, C(f)+E kin (f) < 00} , 

where M > is prescribed. Eventually, we need to restrict ourselves to axially 
symmetric functions in this set: 

J^:={/e^ I f(Ax,Av)=f(x,v), x,veM 3 , 

A any rotation about the x 3 -axisj. 

For the function Q determining the Casimir functional we make the following 
Assumptions on Q: Q, d f Q, d)Q G C([0,oo[xIR), Q>0, Q(0,-)=0 = 
dfQ(0,-), and with constants C, C, F, F'>0, and 0<fc, fc'<3/2, 
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(Ql) Q(f,P)>Cf 1+1 / k , f>F,PeM, 
(Q2) d}Q(f,P)>0, f>0, PelR, 

(Q3) Q(f,-) is increasing on ] — oo,0[ and decreasing on ]0,oo[, />0, 

(Q4) Q(/,0)<C7 1+1/fc ', f<F'- 

The first two assumptions are essential whereas the assumptions (Q3) and 
(Q4) can be modified in various ways, cf. the next section. Note that under 
the assumptions above dfQ(-,P) : [0,oo[— > [0,oo[ is one-to-one and onto for 
any PelR. The steady states which we shall obtain will be of the form 

f Q (x,v) = <l>(E Q -E,P), 

where the particle energy E is defined in terms of the induced potential Uq 
as in 

and E < is some cut-off energy. We should point out that if Q does not 
depend on P which is included in our assumptions — then no symme- 

try assumptions need to be made in the choice of the set Tu nor anywhere 
else. We now state our first main result: 

Theorem 1 For every M > there is a minimizer f G J 7 ^ of the energy- 
Casimir functional Tic- Let Uq denote the potential induced by /q. Then f 
is a function of the particle energy E and angular momentum P as defined 
in /[Ll\) and flTgj, 

f (x,v)=4>(E -E,P) 

where <fi is defined in terms of Q by ( \2.1\ ). The parameter Eq plays the role 
of a Lagrange multiplier and is given by 



Eo:= Jd 1 1 i E + d fQ(fo,P)}fodvdx. 



In particular, f is an axially symmetric steady state of the Vlasov-Poisson 
system with total mass M . 
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The crucial part here is the existence assertion a detailed proof of which 
is given in the appendix. If one computes the Euler-Lagrange condition 
corresponding to our variational problem one obtains the assertions on the 
form of the minimizer. For the details of this tedious but straight forward 
argument we refer to |l(J or [|TT . 



Remark 1. If one is interested only in the existence of a steady state 
the usual approach is to prescribe some suitable function f (x,v) =(p(E — 
E,P), compute the corresponding spatial density po which becomes, via E, 
a functional of the potential Uq, and one is then left with the problem of 
proving that the semilinear Poisson equation 

Af/ = 4vr J '<f>(Eh-±\v\ 2 + U Q ,P(x,vj) dv (2.2) 

has a suitable solution. The non-trivial problem here is to decide which 
solutions have finite total mass and compact support in space. The first 
requirement for </> from the point of view of our theorem is that is a strictly 
increasing function of Eq — E so that a corresponding function Q and the 
Casimir functional can be defined. Of course, for the existence question 
alone this monotonicity condition which in astrophysics textbooks appears 
as a stability condition is not necessary. A mathematical approach which 
makes no such monotonicity assumption can be found in |2~o| . The point is 
that the minimizer which we obtain is not just any old steady state, but a 
nonlinearly stable one. 

Some remarks on stability concepts. It may be useful to briefly 
review various concepts of stability before we state and discuss our stability 



theorem; for a more detailed such discussion we refer to |14]| . Consider a 
nonlinear dynamical system which for short we write as 

x = A(x), 

A being some nonlinear operator on the possibly infinite dimensional state 
space where solutions t^x(t) take their values. Let x be a steady state, 
i. e., A(xq) = 0. 

• The "appropriate" stability approach is certainly one where we can 
deal with our problem directly, i. e., without simplifying or linearizing 
it in any way. The steady state is called nonlinearly stable or Liapunov 
stable if one can guarantee that sufficiently small perturbations of the 
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steady state launch solutions which are arbitrarily small perturbations 
of the steady state for future times. More formally: One must find two 
norms || • ||, ||| • ||| on the state space such that for any (small) number 
e > there always exists another (small) number 5 > such that every 
solution x(t) with \\x(0) — x \\ <5 satisfies \\\x(t) — x \\\ <e for all t>0. 
Some comments are in order: 

— It is desirable but not always possible to use the same norms for 
measuring the initial perturbation and the one at time t; only on 
finite dimensional spaces are all norms equivalent. 

— Sometimes one may not even obtain norms but more general tools 
to measure the deviations with. 

— From a physics point of view it is unsatisfactory if no rule is given 
on how to obtain the S if e is prescribed. From this point of view 
a better result would be an estimate of the form 

|||x(t)-:ro|||<C||a;(0)-a;o||, 

possibly with some explicit constant C. Unfortunately, mathe- 
maticians are not always (not often?) able to provide this. 

— A global existence result which guarantees that solutions exist for 
all time t > at least for initial data close to the steady state is a 
necessary prerequisite and integral part for all the above. 

• Since assessing nonlinear stability in the sense above is difficult, the 
problem is often approached via linearization. By Taylor expansion 
one can compute the linearization of A at x ; formally A\ in — D x A(x ). 
There are now at least two sub-concepts of linearized stability: One is 
to repeat the definition for Liapunov stability given above, but replace 
A by its linearization A\ in . After all, saying that A is nonlinear does 
not forbid A happening to be linear. This is what we want to call 
linearized stability. A somewhat different concept is spectral stability: 
One considers solutions of the linearized equation of the form e xt x and 
tries to find out what the possible eigenvalues A are. If the real parts 
of all these eigenvalues are strictly negative then one calls the steady 
state spectrally stable. Again, some comments are in order: 
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In general there is no guarantee that the linearized problem actu- 
ally has solutions of the form e xt x. If it does, then in order to draw 
any conclusions one must know that there are sufficiently many 
such solutions to get the general solution of the linearized problem 
by superposition (completeness of the eigenfunctions). And even 
if this goes well, there is simply no general result which will allow 
any conclusion on the behavior of the original, nonlinear system, 
unless the system is finite dimensional! 

If one is dealing with a conservative system with some sort of 
Hamiltonian structure (such as the Vlasov-Poisson system), then 
the best to expect as far as spectral stability is concerned may 
be that all the eigenvalues are purely imaginary (due to inherent 
symmetries of the spectrum). In this case no stability follows for 
the nonlinear system, not even in the finite dimensional case. 

Assume one can establish the existence of a growing mode, i. e., 
of a solution of the form e xt x where A has positive real part. Then 
again there is no general result saying that xq is now also nonlin- 
early unstable, unless the system is finite dimensional. However, 
such a growing mode is a valuable first step toward proving a non- 
linear instability result. How difficult it still may be to get to a 
mathematically rigorous nonlinear instability result from there is 
illustrated by ^ [T3 



The upshot from all this is: One should be very careful to draw conclusions 
about stability for nonlinear, infinite dimensional systems from linearization, 
in particular, for Hamiltonian ones. 

All this said let us now return to our minimizer f from Theorem [I]. To 
investigate its dynamical stability we note first that if Jf f — M, 

nc(f)-nc(fo)=d(f,f )-^\\vu f -vu \\l, (2.3) 

where 1 1 - 1 1 a denotes the usual norm on L 2 (IR 3 ) and 

d(f, fo):=ff \Q(f, P) -Q(fo,P) + (E-E )(f- f 



dvdx; 

we are allowed to subtract the term Eo(f — /o) from the integrand since its 
integral vanishes. Since Q(-,P) is convex, the integrand can be estimated 
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from below by 

d f Q(fa,P)+E-E ](f-fa). 

According to Theorem [1], this quantity is zero where fo > 0, while it equals 
(E — E ) f > where f = 0. Thus we see that 



d(/,/o)>o, fe^ 



M- 



Now we note that the left hand side of (|2.3| ) is constant along axially symmet- 
ric solutions. Moreover, and this is the crucial point, the term || V£7 / — VUq ||| 
which seems to have the wrong sign in (|2.3| ) vanishes along minimizing se- 
quences by Theorem AjTJ stated in the appendix. 

Before we state our stability result we point out that if we shift a mini- 
mizer in the x 3 direction we obtain another minimizer. Moreover, we do in 
general not know whether the minimizers are unique up to spatial shifts. As 
discussed in the next section a minimizer must a-posteriori have the following 
additional symmetry property 

f (x 1 ,x 2 ,x* 3 -x 3 ,v 1 ,v 2 ,-v 3 ) = f (x 1 ,x 2 ,x* 3 + x 3 ,v 1 ,v 2 ,v 3 ), x,vE~R 3 , 

for some x^ G IR. If we take without loss of generality x 3 = we refer to this 
as reflexion symmetry. This symmetry propagates along solutions of the time 
dependent problem, and, by restricting ourselves to data with this symmetry, 
we can at least ignore the non-uniqueness due to shifts: 

Theorem 2 For every e > there is a 5 > such that for any solution 1 1— > 
f(t) of the Vlasov-Poisson system with /(0) G C^(IR 6 ) HJ 7 ^ and reflexion 
symmetric, 

d(f(0),fo) + ^\\VU m -VU \\l<8 
implies that for all t>0, 

d(f(t)Jo) + ^\\VU m -VU fo \\ 2 2 <e, 
provided fo is unique or at least isolated in the reflexion symmetric subset of 

•r m- 
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Since after what was said above the proof is simple and instructive we 
include it here: 

Proof. Assume the assertion were false. Then there exist e > 0, t n > 0, and 
/ n (0) G C*(IR 6 ) fl J 7 ^ reflexion symmetric and such that for all nG IN, 

d{f n (0)Jo) + ^ || VC//4Q) - VtfollI < - (2-4) 

but 

/o) + ^ II VU> B w - VC/oIll > 6. (2.5) 

By (0) and 

limH c (fn(0)) = h s M , 

where h s M is the finite, negative minimum of Tic on cf. Lemma A|| Since 
Tic is conserved along classical, axially symmetric solutions as launched by 
/ n (0), since mass is conserved, and since the assumed symmetry propagates, 

YimH c {fn{tn)) = h s M and f n (t n ) G n G M, 

i. e., (fn(t n )) is a reflexion symmetric minimizing sequence for 7i(7 in -^m- 
Up to a subsequence we may therefore assume by Theorem A [I] that 

l|W/„( tn )-V[/ /o ||^0; (2.6) 

note that due to the reflexion symmetry of the minimizing sequence and any 
minimizer the shifts along the X3-axis in the statement of Theorem A[l| must 
be zero in the present situation. It is at this point that the uniqueness or 
isolation of the minimizer f Q is used. Since lim^oo TCc (f n (t n )) = h s M = Hc(fo) 
we conclude by ( p.6|) and (|2.3| ) that 

d(fn(t n )J )->0, rwoo, 

and we arrive at a contradiction to ( |2.5|) . □ 

Remark 2. If the minimizer is not isolated, i. e., if arbitrarily close 
to it with respect to the "distance" used in the theorem there are other 
minimizers which do not result from a simple translation in space, then we 
obtain a stability result of the following form: 
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Let A4m C J-'m denote the set of all minimizers of Tic i n J~m- Then for 
every e > there is a 5 > such that for any solution ti—>f(t) of the Vlasov- 
Poisson system with /(0) G C^H 6 ) njFf^ , 



inf 

fo&M M 



^(/(0),/o) + ^||V[/ /( o)-VC/o||^ 



implies that for all £>0, 



<5 



inf 

h&M M 



d(f(t),f ) + — \\VU m -VU h 



<e. 



Remark 3. Although we only showed that d(f,f )>0 for /GJF M , one 
may think of this term as a weighted L 2 -difference of / and / , and by a 
Taylor expansion a stronger and more explicit estimate can be obtained if 
djQ is bounded away from zero. 

Remark 4. The restriction /(0)GJF M for the perturbed initial data 
would be acceptable from a physics point of view: A perturbation of a given 
galaxy, say by the gravitational pull of some outside object, results in a 
perturbed state which is just a rearrangement of the original state, in par- 
ticular, its mass remains unchanged. But in general symmetry properties 
of the steady state will be destroyed by such a perturbation. However, we 
need that /(0) is axially symmetric so that C is conserved along the resulting 
axially symmetric solution f(t). So while this symmetry restriction for the 
perturbation seems necessary within the present mathematical framework it 
is quite undesirable from a physics point of view. If the steady state does 
not depend on P and thus is isotropic, then no symmetry restrictions are 
necessary anywhere in our results, and the stability is with respect to quite 
general perturbations. 

Remark 5. Our stability result is of the nonlinear Liapunov type dis- 
cussed above in general, but it suffers from the defect that given e it does 
not say how 5 needs to be chosen. The reason for this is that our proof is by 
contradiction and relies on a compactness result. 



3 Examples and further properties 

In this section we investigate some additional properties of the resulting 
steady states, keeping the discussion somewhat informal. First we present 
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some examples for functions Q which satisfy our assumptions or possible 
variations of them: 

Examples. A simple class of examples for functions Q which satisfy the 
assumptions (Q1)-(Q4) is given by 

Q(f,P) = f 1+1/k g(P), />0, PeM, (3.1) 

where < k < 3/2 and g is positive, continuous, bounded, bounded away from 
zero, increasing on ] — oo,0[, and decreasing on ]0,oo[. Note that in this case, 

f (x,v) = C(E -E) k g(P)- k 

on its support. Examining the proof of Lemma A [3] in the appendix, which 
is the only place where the assumptions (Q3) and (Q4) enter, one can see 
that these assumptions can for example be replaced by 

(Q3') Q(Xf 1 P)>X 1+1 / k 'Q(f,P), f>0, 0<A<1, PeJR, with 1/2 < k' < 3/2, 

(Q4') Q(f,P)<C'f 1+1 / k ", f<F', \P\>P , with constants C>0, F'> 
0, P >0, and 0<k"<3/2. 

Examples which satisfy the assumptions (Ql), (Q2), (Q3'), (Q4') but not 



the original ones are provided by (|3.1| ) with l/2</c<3/2 and g which has 



all the properties stated above except the monotonicity. 

Regularity and boundary condition. Since by Theorem [l] f is a function 
of the quantities E and P which are constant along particle trajectories we 
are justified to call fo a steady state provided Uq is sufficiently regular to 
allow for the definition of particle trajectories to begin with. Now by con- 
struction AUo = 4:Tcp on IR 3 , at least in the sense of distributions, and from 
the very construction of fo certain integrability properties of po follow, cf. the 
appendix. Applying the usual Sobolev space arguments one can eventually 
conclude that Uq is twice continuously differentiable and lim| x .|^ 00 f/o(x) = 0, 
cf. [O, Thm. 3] for the technical details. 



Finite mass and compact support. By construction, the steady states have 
finite mass M>0 which we prescribe by our mass constraint. To continue 
we note that by Theorem |1] and a change of variables, 



po(x) = J f (x,v)dv 

= 2ti / t <p{E -E,r(x)p)dpdE (3.2) 

JUo(x)J-y/2(E-Uo(x)) 
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where Uq(x) < E , and po is zero else; here r{x) :— yxf + a;|. Assume that 
Hmui^ooC/o^) = 0, see the discussion above. Then clearly the cut-off energy 
cannot be positive, i. e., E <0, since otherwise we get infinite mass "at 
spatial infinity" . Since Uq < everywhere, (|3.2|) shows that p wn l have 
compact support if and only if E <0. An additional, sufficient condition 
which implies this is 

fd f Q(f,P)<3Q(f,P), f>0, PeTR, (3.3) 

which holds for example for Q(f,P) = f 1+1 / k g(P) with l/2<£;<3/2 and 
some function g. To see this we rewrite the formula for Eq from Theorem [1} 



1 

Eo= M 



d 

dl. 



E kin (f )+2E pot (fo) + J J f d f Q(fo,P)dvdx . 
Now we note that for solutions of the Vlasov-Poisson system, 

■ JJx-vf = 2E Un (f) + E pot (f) 

so that for a steady state the right hand side is zero. Thus by (|3.3|) , 

E k Ufo)+2E pot (fo) + J J f d f Q(fo,P)dvdx < 
E kin (/ )+2E pot (/ )+3 J J Q(f ,P)dvdx = 

3E kin U ) + 3E vot (f )+3C(f ) = 3W G (/o) < 

by Lemma A |^ (a). An alternative condition which also guarantees compact 
support is the following: 



(f>(E -E,P)<Cx(E -E) kl , E^-oo, PeTR, 
0(E o -E,P)>C 2 (E o -E) k \ E^E -, PEM 



(3.4) 

for positive constants Cj, C 2 and < ki,k 2 < 3/2. This assumption also im- 
plies that £ <0, cf. @, Thm. 3]. 

In the spherically symmetric case our approach provides an explicit bound 
on the radius of the steady state, cf. 0. Note that in this case one can under 
appropriate assumptions decide whether Uq crosses the cut-off energy level 
E by direct examination of the semilinear Poisson equation (|2.2|), since this 
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equation reduces to an ordinary differential equation with respect to the 



radial variable, cf. p5| . However, in the axially symmetric case this is no 
longer true, and the corresponding analysis of the genuine partial differential 
equation (|2.2|) would be much more difficult. 

Symmetry. By construction the steady states which we obtain are axi- 
ally symmetric, but so are the spherically symmetric ones whose existence 
was known already In order to make sure that we have found qualitatively 
new steady states we must show that they are in general not spherically 
symmetric. To see this we take Q such that <P{Eq — E,-) is not constant on 
any neighborhood of P = 0. We claim that neither po nor Uq are spherically 
symmetric in this case. Indeed, if po were spherically symmetric the same 
were true for Uq. If Uq were spherically symmetric we take some 2G H 3 with 



r{x) = Jx\ + x%^0 small and take AeSO(3) such that Ax— (0,0, \x\), hence 
r(Ax) = 7^ r(x) but Uq{Ax) = Uq(x). Inserting this into the formula ( |3.2| ) for 
Po we see that in general p (Ax) ^ p (x) so p is n °t spherically symmetric. 

Indeed, this is a bit more than just saying that fo is not spherically 
symmetric. To see this, take a spherically symmetric steady state fo with 
induced density po, potential Uo, and particle energy E— \v \ 2 /2 + Uq(x). Let 
g (x,v) =ip(E,P) be any function which is odd in P, and such that fo + go > 0. 
Since P is not invariant under general rotations neither is fo + go, but the 
fact that if) is odd in P implies that pf 0+go = p and the same holds for the 
potential. So this trivial construction gives an axially symmetric steady 
state where the phase space density is not spherically symmetric but the 
macroscopic quantities po and Uq are. Our steady states are in general not 
of this trivial type. 

Another symmetry issue refers to the dependence of the minimizer on the 
variable 23. Given a minimizer f of which we do not yet know any sym- 
metry with respect to 23 we can do a symmetric decreasing rearrangement 
of fo with respect to 23 while keeping all other variables fixed. Denote this 
rearrangement by /q*; as to the rearrangement concept we refer the reader 
to 116, Ch. 3]. Obviously, the rearrangement does not change the kinetic 
energy, and neither does it change the Casimir functional, since P does not 
depend on 23. By [0, Thms. 3.7, 3.9] it can at most decrease the potential 
energy, with equality iff fo is already symmetric and decreasing in 23 up to 
a shift in 23. But since fo minimizes the energy-Casimir functional it must 
posses this symmetry, and without loss of generality we can assume that fo 
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is reflexion symmetric in X3. The corresponding symmetry in the variable V3 
follows from the form which /o must have due to Theorem [l]. 

Stationary versus static solutions. If instead of the Vlasov-Poisson system 
we consider a self-gravitating fluid as described by the Euler-Poisson system 
then every so-called static solution, i. e., every steady state with vanishing 
velocity field, is spherically symmetric, cf. ||15|| . This turns out to be false 
for the steady states of the Vlasov-Poisson system which we obtain. By 
definition the velocity field equals jo/po on the support of po where the mass 
current density jo is given by 

r rE ry/2(E-U (x)) 

jo(x) = / vf (x,v)dv = 2ir / / s<j)(E - E,r(x)p)dpdE e t (x) 

J JUo(x)J-y/2{E-Uo(x)) 

with e t (x) := (— X2,xi,0)/r(x); on the X3-axis the velocity field vanishes. Ob- 
viously, if is even in P then jo vanishes identically, so there exist static 
solutions which are not spherically symmetric among the ones we obtain. 
On the other hand, if, say, <j>(E,P)><f>(E,-P) for all E>0 and P>0 the 
velocity field is non-trivial and corresponds to an average rotation about the 
axis of symmetry in the counterclockwise direction. 

Less regular dependence on P. When examining the appendix the reader 
may notice that at no place we really make use of the fact that all the 
functions under consideration can be taken to be axially symmetric. This 
is because as opposed to spherical symmetry, exploiting axial symmetry for 
example in the estimates for the potential is technically quite unpleasant. 
However, it is possible and even necessary if one wishes to study examples 
of the form 

Q(f,P) = f+Vk p -2i/k^ / > 0, P e IR (3.5) 

which for / 7^ do not satisfy the assumptions on Q which we stated above. 
The corresponding steady states would take the form 

fo{x 1 v) = C{E -E) k P 21 

which is analogous to the classical, spherically symmetric polytropes, except 
that the square of the third component of the angular momentum replaces 
the square of the modulus of the angular momentum. Clearly, ( |3.5| ) satisfies 
the crucial convexity condition (Q2), provided k>0. When we examine the 
scaling arguments in Lemma A [5] we find that they go through, provided 

3 

1>-1, 0<k<l + -. 

2 
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The place where significant changes become necessary are the estimates in 
Lemma A |l| and Lemma A || Clearly, one can no longer control ||/||i+i/k 
and ||p/||i+i/ n in terms of E kin (f), C(f), and M. Instead, one gains control 
of the weighted norms 

J J p- l / k (x,v)f 1+1 / k (x,v)dvdx, J r- 2l l n (x)p) +1/n {x)dx 

where r(x) = \jx\-\-x\. We conjecture that the weight factor can be dealt 
with in the estimates for the potential energy etc. if one makes proper use of 
the axial symmetry of the potential. 

4 Appendix 

We now give the proof of the existence part in Theorem [I] and the additional 
results used in the proof of Theorem |2|. First we collect some estimates for 
Pf and Uf induced by an element /£ JF M . These estimates make no use of 
symmetry and rely on the assumption (Ql). Their main point is to see that 
Tic is bounded from below on Tu and to establish certain bounds along 
minimizing sequences of Tic- Constants denoted by C are positive, may 
depend only on Q and M, and their value may change from line to line. By 
II • |L we denote the usual L p -norm of functions over H 3 or 1R 6 as the case 
may be. 

Lemma A 1 Let n := k + 3/2 < 3. Then for any f G Tm the following holds: 

(a) J J f 1+1/k (x,v)dvdx<C{l + C{f)) 

(b) / p) +l/n {x)dx<C (l + £U/) + C(/)) 

(c) U f EL 6 (JR 3 ) with VU f EL 2 (JR 3 ), the two forms of E pot {f) stated in 
Section 2 are equal, and 

J \VU f \ 2 dx < c \\ Pf \\l /b < C (1 + E kin (f) +c(f)) n/3 . 

Proof. Part (a) is a direct consequence of (Ql), splitting the v-integral in 
the definition of pf into small and large v's and optimizing the split yields 
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(b), while the extended Young's inequality and interpolation together with 



(b) imply (c). For details cf. [ 10 1 or [11]. □ 



The estimates above have an immediate but important consequence: 

Lemma A 2 The energy- Casimir functional Tic is bounded from below on 
T M , i- e., 

h M := inf H c (f) > -oo, 

and along any minimizing sequence of Tic in Tm the quantities E^i Q (f), C(f), 
H/Hi+i/fc, and ||p/||i+i/n are bounded. 

Proof. Lemma A [1| yields the estimate 

H c (f)>EUf)+C(f)-C(l + E kin (f)+C(f)) n/ \ feT M , 

and since n < 3 the assertions follow. □ 

Lemma A [l] and Lemma A ^| remain valid if we replace Tm by its axially 
symmetric subset T M . Of course the infimum of Tic on this smaller set may 
be larger; we denote 

h S M := MHc(f). 

The assumptions (Q3) and (Q4) determine the behavior of He under 
scaling transformations which we use to show that hu is negative and to 
relate the /ijw's for different values of M: 

Lemma A 3 (a) LetM>0. Then — oo < h M < 0. 

(b) ForO<M<M, 

h M >(M/Mf /3 h M . 

Proof. Given any function /, we define a rescaled function f(x,v) = 
af(bx,cv), where a, b, c> 0. Then 

J J fdvdx = a(bc)~ 3 J J fdvdx, (4.1) 
i. e. / G Tm iff f^T M where M = a(bc)~ 3 M. Next 

Hc{J) = ab- 3 c- 5 E kin (f)+a 2 b- 5 c- 6 E pot (f) 

+ (bc)- 3 [ [ Q{af{x,v),{bc)- 1 P(x,v))dvdx. (4.2) 



18 



To prove (a) we fix a function f £T\ with f <F' and let a = M(bc) 3 so that 
/ G Tm- Then by (Q3) and (Q4) and with positive constants Ci,C 2 ,C 3 which 
depend on /, 

Hc(f) < Mc- 2 E k in(/)+M 2 6 J E; pot (/) + Ma- 1 J J Q(af(x,v),0)dvdx 

<c lC - 2 -c 2 b+c^ k \ 

provided a< 1 so that (Q4) can be applied; note that E pot (f) <0. We want 
the negative term to dominate as 6—^0, so we let c = b~ 1 ^ 2 for some 7>0. 
Then a = Mb 3(l ^l 2 \ and 

Uc{f) < C X W - C 2 b + C 3 M 1/k 'b 3(1 -~ t/2)/k ' . 

Since 0</c'<3/2 we can fix 7G]1,2[ such that 3(1 -y/2)/kf> 1. For b> 
sufficiently small Ti.c(f) will then be negative and a = M6 3 ^ 1-7 / 2 ^ < 1 as 
required. This proves part (a) of the lemma. To prove (b) we choose a = c = 1 
and b = (M/M)" 1//3 so that the mapping Tm Fm-i f^f^ s one-to-one and 
onto and 6 1 < 1. By (§3), 

nc(f)=b- 3 E kin (f)+b~ 5 E P M)+b- 3 JjQ(f(x,v),b- 1 P(x,v))dvdx 

> b- 5 E kin (f)+b- 5 E pot (f)+b- 5 J J Q(f(x,v),P(x,v))dvdx 
= b- 5 H c (f); 

we multiplied the two positive terms by b~ 2 < 1 and used the monotonicity 
of Q in the P- variable which we required in (Q3). By the choice of b and the 
definitions of hu and h M the proof is complete. □ 

It is again obvious that the assertions of the lemma remain valid if we 
restrict ourselves to the axially symmetric functions in Tm- The reader might 
also check that the scaling arguments work under the assumptions (Q3') and 
(Q4') from Section 3 as well. 

Next we provide a splitting estimate which will be used to show that 
along a minimizing sequence the mass cannot escape to infinity; here and in 
the following we denote for < R < S < oo, 

B R := {xGlR 3 ||x|<ii!}, 
B R>S ■= {sGlR 3 |i?<|x|< 1 S}. 
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Lemma A 4 Let f&J~M- Then 

The proof follows from splitting the potential energy into three parts ac- 
cording to \x — y\ < 1/R, l/R< \x — y\ <R, and R > \x — y\; for the details we 
refer to (Tl| or |23|, Lemma 3] The splitting estimate above has the following 



important consequence for minimizing sequences: 

Lemma A 5 Let (/i)cJF M be a minimizing sequence of Tic- Then there 
exist <5o>0, i?o>0, io^lN, and a sequence of shift vectors (a«)cIR 3 such 
that 

\ fi(x,v)dvdx>5 , i>io, R>Rq. 

Jai+B R J 

If (fi)cJ-'M, i. e., the minimizing sequence consists of axially symmetric 
functions then one can choose ai = (0,0, Zi) with Zi G IR, i. e., only shifts along 
the axis of symmetry need to be admitted. 

Proof. By Lemma A 0, (||p/, ||i+i/n) is bounded. By Lemma A |3] (a) we 
have 

E pot (fi)<Hc{fi)<^h M <0, i>i Q , 

for a suitable i G IN. Thus by Lemma A [| there exist S > 0, R > 0, and a 
sequence of shift vectors (a^) cH 3 as required. 

Now assume that the functions fi are axially symmetric, and let = 
(di,Zi) with a,i G IR 2 and G IR. Suppose we can show that the sequence (ai) 
is bounded. Then we can pick some R > such that 

(d h 0) + B Ro cB R , ieIN 

and thus 

a* + 5^(0,0,^) + ^, zgIN, 
which implies our assertion in the symmetric case: For R>R, 

fi(x,v)dvdx> / / fi(x,v)dvdx>5o, i>io- 

l(0,0,Zi)+B R J Jai+B Ro J 
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So it suffices to show that (cij) is bounded. Assume the contrary and choose 
some integer AG IN such that N5q>M. If \a~i\ is sufficiently large then 
simple geometry tells us that there exist rotations Ai,...,A N about the x 3 - 
axis such that the balls A/ C a i + -Br , k = l,...,N, are pairwise disjoint. By 
axial symmetry of 

/ fi(x,v)dvdx= / fi(x,v)dvdx>5o, k = l,...,N, 

JA k ai+B Ro J Jai+B Ro J 

and since the balls on the left hand side are pairwise disjoint, 

N 

/ / fi(x,v)dvdx>^2 / / fi(x,v)dvdx> N 5 > M 

J J k=l JA k a i +B Ro J 

which violates the mass constraint and gives the desired contradiction. □ 

We will also need to exploit the well known compactness properties of the 
solution operator of the Poisson equation: 

Lemma A 6 Let (pA C L 1+1 / n (]R 3 ) be bounded and 

Pi^Po weakly in L 1+1/n (lR 3 ). 

(a) For any R>0, 

VU 1brPi ^VU 1brP0 strongly m L 2 (IR 3 ). 

(b) If in addition (pi) is bounded in L 1 (IR 3 ) ; poE^fJR 3 ), and for any e>0 
there exists R>0 and io G IN such that 

\pi(x)\dx <e, i>io 

x\>R 

then 

VU Pi ^VU pa strongly in L 2 (1R 3 ). 

For the details of the proof we refer to |TTJ or |23], Lemma 5]; it relies on 
the fact that over bounded sets the solution operator to the Poisson equation 
is compact in the appropriate spaces and the fact that in the context of the 
lemma we can make the contribution to the potential energy coming from 
outside sufficiently large balls small. 

We are now ready to prove the existence of a minimizer, more precisely 
we prove: 
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Theorem A 1 Let (/;) C be a minimizing sequence of Tic- Then there 
exists a sequence of shifts (zi) C IR and a subsequence, again denoted by (ft), 
such that for any e > there exists R>0 with 

[ fi(x,v)dvdx>M-e, i G IN, 

J(0,0,Zi)+B R 

T£ :=£(• + (0,0,*), 'Wo weakly in L 1+1 / k (WL 6 ), woo, 

and 

/ fo>M-e. 

JB R 

Finally, 

Vt/r/w Vt/ strongly in L 2 (IR 3 ), woo, 
and fo EJ 7 ^ ^ s a rninimizer of Tic- 

Proof. We split / G Tf^ into three different parts: 

/ = 1 B Rl xtf/ + 1 B Rl;R2 xtf/ + 1 B R2 , 00 xK 3 / = : fl + h + /s5 

the parameters Ri < R 2 of the split are yet to be determined. The induced 
spacial densities are denoted by p\~, their masses by and the induced 
potentials by £4, k — 1,2,3. With 



// m := / / — | 1 — , /,m = l,2,3, 

we have 

HcU) = Wc(fi) + ndh) +n c {h) - in - hs - hs- 

If we choose R 2 > 2R X then 

3 ~ R 2 

Next, we use the Cauchy-Schwarz inequality, the extended Young's inequal- 
ity, and interpolation to get 



1 

I\2 + -^23 — ~r~ 

4:71 



Jw(U 1 + U 3 )-WU 2 dx <C||pi+p 3 ||6/5||Vi7 2 || 2 

<^llp|lS n + vi /6 llvc/ 2 || 2 . 
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Using the estimates above and Lemma A |3| (b) we find 



M J \MJ \M 
+G(i? 2 - 1 +||p||K /6 ||V[/ 2 || 2 

c 

-l_l_ll„l|( n + 1 )/6| 



< — (MiM 2 + MxM 3 + M 2 M 3 ) /if, 



+c(i? 2 - i +iipii^;rnvf/ 2 ii 2 

< Ch 9 M M x M^+C (r, 1 + \\ P \\tV/n 6 l|V£/ 2 ||a) ! (4.3) 

observe that by Lemma A ^ (a) /if f < and that constants denoted by C 
are positive and depend on M and Q, but not on Ri or R 2 . We want to 
use ( |Q| ) to show that up to a subsequence and a shift M 3 becomes small 
along any minimizing sequence for i large provided the splitting parameters 
are suitably chosen. 

So let (fi)cJ-M be a minimizing sequence and define Ro, 5q, and (zj) 
according to Lemma A [5]. The shifted sequence (Tfi) is again minimizing so 
by Lemma A Q we can pick a subsequence, again denoted by (/j), with a weak 
limit as stated in the theorem. Now choose Ri > R Q so that by Lemma A [5], 
M hl >5 for % large. By©, 

(4.4) 

where U^i is the potential induced by f it i which in turn has mass Mij, ie 
1NU{0}, and the index Z = 1,2,3 refers to the splitting. Given any e>0 we 
increase Ri > Rq such that the second term on the right hand side of (|4.4j) 
is small, say less than e/4. Next choose i? 2 >2R\ such that the first term is 
small. Now that R\ and i? 2 are fixed, the third term in ( |4.4| ) converges to 
zero by Lemma A || (a). Since (Tfi) is minimizing the remainder in (|4.4| ) 
follows suit. Therefore, for i sufficiently large, 

/ [Tf^M-M^M-i-ChffSoyh. (4.5) 

J {0,0,Zi)+B R2 J 

Clearly, fo>0 a. e., and fo is axially symmetric. By weak convergence we 
have that for any e > there exists R > such that 



-Chf 4 6 Q M it3 <— + C\\VU 4 2 + C\\VU h2 -VU 0i2 \\ 2 + n c (Tf l )-h s M 



M> [ { f dvdx>M-e 
Jb r J 
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which in particular implies that foE^fJR 6 ) with / fodvdx = M. The func- 
tional C is convex, so by Mazur's Lemma and Fatou's Lemma 

C(/ )<limsupC(T/-). 

i— >oo 

The strong convergence of the gravitational fields now follows by Lemma A |6] 
(b), and in particular, 



n c (fo)<^sn P Hc(Tf i ) = h 



M 



so that /o G Tfa is a minimizer of He- n 
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